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The propagation of low-momentum baryons through QCD matter at low isospin density is studied 
using chiral perturbation theory. When the isospin chemical potential exceeds a critical value, 
the dispersion relation for the lowest nucleon mode becomes anomalous near zero momentum — 
increasing momentum yields decreasing energy — so that the momentum for the lowest energy state, 
Pmin, becomes nonzero. This can be interpreted as a new phase of QCD, with p m in serving as 
the order parameter and the kinetic mass (defined as the inverse of the second derivative of the 
energy with respect to the momentum at p m in) serving as the susceptibility. A lowest-order chiral- 
perturbation-theory calculation yields a critical isospin chemical potential of ~ 285 MeV. Since this 
is small compared to the chiral-symmetry-breaking scale, corrections to it are likely to be modest. 



Recently, Son and Stephanov have considered QCD in a phase with nonzero isospin density fit] . Such a system is 
of interest for a number of reasons. It is an intrinsically nonperturbative system which is sensitive to much of the 
complex dynamics of QCD. In the limits of both very low and very high isospin density, there are analytical techniques 
which, given our present state of understanding of the theory, should be both valid and computationally tractable. 
At the same time the functional determinant for the system in Euclidean space is manifestly real and positive [0. 
Accordingly, one can numerically simulate QCD at nonvanishing isospin density, p f , (or more precisely nonvanishing 
isospin chemical potential, fir) using standard lattice Monte Carlo techniques, unlike the case of QCD at finite baryon 
density j^]. Thus, even if it is hard to envision a situation in which one could study this regime experimentally, 
the prospect of lattice results means that it may serve as a useful testing ground for various theoretical approaches 
to QCD. The system may also provide insights into some of the properties of neutron matter — a system with both 
non-zero isospin density and baryon density. Moreover, the phase structure of the theory is interesting in its own 
t-H ■ right. At low, but non-zero, isospin density the medium is a superfluid — a pion condensate. The phase structure is, 
in fact, more rich. For example, if one generalizes isospin to three flavors as was recently considered by Kogut and 
Toublan one finds a kaon condensed phase which is separated from the pion condensed phase by a first order phase 
transition. In this paper we show that even if one restricts one's attention to two flavors and relatively low isospin 
densities, the phase structure of the theory is quite nontrivial — in addition to the ordinary and pion-condensed phases 
discussed in ref . , a new phase with radically different nucleon properties emerges at a critical isospin density. 

Here we will focus on the low p t regime. The theoretical tool used by Son and Stephanov to study this regime is 
quite conventional — chiral perturbation theory (%PT) The essential point is that an isospin chemical potential 
couples to the quarks in QCD like a time-like photon. Thus, by simply imposing gauge invariance at the level of the 
effective theory one has included the leading effects of the chemical potential. Clearly, the isospin chemical potential 
will have to be restricted to values that are small compared to the characteristic hadronic scale for xPT to be valid. 
As noted in ref. if the isospin chemical potential is larger than m w one is in a chirally nontrivial phase. The 
essential physical point is quite simple. For \i t > it becomes energetically favorable to condense pions out of the 
vacuum. Pions continue to condense until the pion-pion interactions (which can be described accurately by xPT) 
yields a repulsion equal to the attraction due to the chemical potential. Thus as the isospin chemical potential is 
increased from zero, the isospin density remains zero until p, = m n at which point p 1 increases smoothly from zero. 
From the perspective of xPT, which is based on a nonlinear sigma model, the pion condensation described as a chiral 
rotation from the isoscalar (a) direction toward one of the pion directions. 

At first sight it might appear that this pion-condensation transition is the only one which can be described in 
a regime where xPT is valid. After all, it is necessary to have an order parameter which distinguishes the pion 
condensate from any putative new phase. In the context of the chiral Lagrangian — a simple theory of interacting 
pions — there is no obvious candidate for such an order parameter aside from the chiral field itself. However, pions are 
not the only possible hadronic probes of this state of matter: one can send nucleons through this matter to probe its 
structure. 

The nucleon dispersion relation can be used to define a new phase. For nucleons in free space, Lorentz invariance 
completely determines the dispersion relation: E 2 — p 2 + Al 2 . However, an isospin chemical potential breaks Lorentz 
invariance since it couples to the zeroth component of a four-vector. Thus, the nucleon dispersion relation in this 
medium can be nontrivial. At sufficiently low fi^ both above and below the pion condensation transition, the 
nucleons have a normal dispersion relation — the energy increases with increasing momentum. However, when |/i f | is 
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above a certain critical value, ^ cr , the situation is fundamentally different; for low values of the nucleon momenta the 
dispersion relation for the lowest lying nucleon excitation is anomalous in the sense of anomalous dispersion in optics 
— increasing the momentum leads to lower energy. This defines our new phase. It can be described by a simple 
order parameter, 

© = pZST > (!) 

where P^^ on is the magnitude of the momentum of the lowest-lying excitation of nonzero baryon number. For 
I/it/ 1 < /Li cr , O = while for |/i/| > /x cr , 6^0; thus O serves as an order parameter in the traditional sense of 
distinguishing between two phases. As will be shown below, the critical chemical potential, /i cr , can be estimated 
using xPT. At lowest order in the chiral expansion the critical chemical potential satisfies 

2M N mln CI = l 



(1 - g 2 ) m\ + g 2 

which yields an estimate of fi CI ~ 285 MeV for the transition. Chiral corrections to this are of order Q 2 /A 2 , so that 
one might expect corrections to alter this value by ~ 20%. 

We begin with a brief recapitulation of the arguments of rcf. jl| for the low \i t region. As all dynamical scales in the 
problem are much less than typical hadronic scales one can describe the system by an effective theory for interactions 
among pions, the only light degrees of freedom in the problem. The Lagrangian for the purely pionic part of the 
theory is simply that of the nonlinear sigma model. At leading order it is 

C = \fl Tr [V M C/V^t + 2ml Re U] , (3) 

where U 6 SU(2) and can be parametrized by U = cxp fiy^V The theory at this order contains two parameters: 

the pion mass, w 139 MeV, and the pion decay constant /„. w 93 MeV; the covariant derivative V incorporates 
couplings to gauge fields. 

At the QCD level the isospin-chemical-potential term is simply \niV IM q'y ,J 'Tzq, where is a four velocity that 
describes the motion of the resulting medium; in its rest frame = (1,0,0,0). Thus, the isospin chemical potential 
couples directly to the conserved vector-isovector current, and at the level of the effective theory, the coupling of the 
isospin chemical potential must be the same as the coupling of a U(l) isovector gauge field. 0. Hence we use 

V ll U = d li U + i[U,V li ], (4) 

with Vfj, — ^HjV^ts = ^HjT^Sq^. Using this lowest order Lagrangian one expects predictions of pionic observables and 
thermodynamic quantities to be accurate up to corrections of relative order Q 2 /A 2 where Q ~ m^,///, p (where p is 
the external momentum of a probe and A ~ 1 GeV is a characteristic hadronic scale). 

The next step is to determine the U which minimizes the energy. The ground state is translationally invariant and 
static so one can simply use a constant U in the Lagrangian of eq. (|^) and equate it to the negative of an effective 
potential, 

f2 2 fl 2 

V eS (U) = I^Tt(t 3 Ut 3 U^ - 1) - ^|^ReTr (£/). (5) 

It is clear that first term in eq. (|J) lowers the energy if U aligns along the t\ or directions while the last term 
favors the isoscalar direction. Accordingly, the minima can be found using an ansatz which allows for general rotation 
between the isoscalar and t± or T2 directions, 

U = cos a + i(n cos cj> + t 2 sin 4>) sin a , (6) 

for which the effective potential takes the form 

Veff(a) = (cos 2a - 1) - f 2 mucosa . (7) 

Note that V e s depends on a but not on <f>. Minimizing with respect to a one sees that for \fj, x \ < , the minimum 
occurs at a = (U = 1) so that one is in the normal phase. For |/ij| > m^, the minimum occurs at 0] 

9 

Tfl 

cos a = — £ . (8) 
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FIG. 1. The nucleon dispersion relation for values of \ii below and above [i< 



which implies that the state has chirally rotated to a new pion-condensed state. 

Now let us turn to the main subject of the present work, the propagation of nucleons through a pion-condensed 
phase. The interaction between nucleons and the chiral field can also be described using X-PT. The lowest order 
Lagrangian describing the interaction of nucleons with pions in x?T is given by Q 

C^n = * (ip- M + iff^Ts^) * , (9) 

where u 2 — U and 

= + ^[u\ d M u] - ^i(u^V M u + uV^v)) , 
Uf_ L = i(vJV - ■uV AI u t ) , 
S7 fJ ,u = d ll u + i[u,Vf l ]. (10) 

We note that the interaction between nucleons and pions considered by Son and Stephanov [Q was based on a nucleon 
field that transforms linearly under chiral transformations rather than one that transforms nonlinearly as considered 
here. In fact, it is straightforward to show that the two are equivalent — that is, they are related to each other by a 
change of variables — provided that one takes g A = 1. The form used in eq. (^|), the standard form used in baryon 
%PT, has the advantage of being able to incorporate g A ^ 1 in a straightforward way. As detailed above, the isospin 
chemical potential acts like a vector field: V" M — ^M t 3^m0- Inserting this into the Lagrangian and varying with respect 
to yields the following Dirac equation: 

[i0 - M + i/i 7o +u^t 3 u) + g Alb \{uT 3 v) -v)t z u))\ * = 0. (11) 

With the ansatz of eq. (||) for the pion field, and a plane-wave state for the Dirac field, this gives 

\E — a ■ p — M70 + \\i (cosqt 3 + 5^75 sin a (cos0t 2 — sin^ri))] = . (12) 

The solutions are eigenstates of helicity, and there are four doubly-degenerate solutions for each value of the momen- 
tum. The two positive-energy solutions are given by 




Although the pion condensate breaks parity, it also breaks CP and thus conserves charge conjugation. Thus, the 
energies of the positive- and negative-energy solutions are equal up to an overall minus sign. The pion-condensed 
phase breaks ^3 since the condensate is aligned in one direction in the t\ — r 2 plane. Thus, in this phase the solutions 
of the Dirac equation do not correspond to states of well-defined isospin — rather they are admixtures of -Z3 = ±1/2. 
Only for a = are the solutions eigenstates of I3 . 

Ref. |ij considered the "nucleon mass", defined as the minimum energy needed to produce a baryon- number-one 
excitation, implicitly taken to be at p = 0. Their expression corresponds to the first two terms in the expansion of 
eq. ( |l3| ) in powers of 1/M. However looking at the full dispersion relation, we see that for a close to ir/2 the point 
p = is no longer the minimum for the lower branch. Expressing a in terms of [i 1 through eq. (^), we find a critical 
value of already given in eq. (0), above which the minimum energy is at a non-zero value of p. Fig. |l| illustrates 
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FIG. 2. The value of p which minimizes the energy as a function of fij, eq. (|l4|). 



the qualitative difference between the two regions. We plot E± for two choices of fXj, one below the critical value and 
one above. 

As noted earlier, p m ; n functions as an order parameter. For > /i cr it is given by 



1 / (g|Mf + (l-gl)mjy - ANPmjp 2 

2 V m? (alf4 + (i - 9a) m i) 



(14) 



This is plotted in fig. g. 

An alternative definition for the mass to that considered in ref. jj]] is the kinetic mass, which for any mode is defined 
by 



Jw kin — 



d 2 E 



dp 2 



(15) 



Clearly for the free-space dispersion relation the mass defined by the minimum-energy excitation and the kinetic mass 
are equal, but in medium this need not be so. The kinetic mass is important in that it is the inertial parameter which 
enters in an effective Schrodinger equation describing the lowest-energy excitations of this system. The kinetic mass 
for the lower of the two nucleon modes is 



Mml ^/(l- 5 2) m 4 + g 2^i _ ^Mmlfij + 4M 2 /i? 



kin 



2Mm2/i, - (l-g\)m* 

g A M((l-gi)mt+girf 



9>i 



\m < Mcr , 
\[J<l\ > Mcr ■ 



(16) 



It is plotted in fig. [| One striking feature is that Mki n develops a pole, diverging at the critical value fi CI . This critical 
value occurs when the denominators of eq. ( |l6| ) vanishes. 

Figures ^| and |^ look very much like a standard second order phase transition, with the kinetic mass diverging at the 
transition in the manner of a typical susceptibility and p m in rising from zero as a fractional power in the manner of 
a typical order parameter. In a way, however, this phase transition is rather unusual. Generally at phase transitions 
there are discontinuities in quantities associated with bulk energetics such as specific heats. In the present case we 
have no evidence of such effects. The basic reason is very clear: the discontinuities associated with the phase transition 
reported here all involve baryon properties. To the order at which we have calculated, the pionic degrees of freedom 
are unaffected by the transition and the pions dominate the energetics. To the extent that there are discontinuities 
in quantities associated with energies they must involve changes in possible nucleon-antinucleon loop effects which 
might altered discontinuously when the nucleon dispersion relation is altered. However, any such possible effects are 
too small to be calculated in %PT. 

We note that the value of the critical chemical potential from Eq. 2 is of order (M N m 2 T ) 1 /' i , and numerically equal to 
285 MeV~ 2771^. Thus in practice /i cr is comparable in size to and well below M N , so we are in the usual regime of 
applicability of xPT. The leading corrections to our result should come from the inclusion of the higher-order terms in 
the pion- nucleon Lagrangian which are parameterised by the low-energy constants Q. In fact to order l/M N the sole 
effect of these is to add a term (4cim 2 — (ca + c^)^, 2 ) sin 2 a to the nucleon energy. Not only is this small (58 ± 1 MeV 
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at p = p CI using the second-order determinations of the Cj Q), it is also momentum-independent. Thus it cannot 
affect /i cr at all, and corrections start only at third order. 

Finally, we turn to a question of considerable significance in nuclear physics, namely neutron matter, which is QCD 
at finite baryon chemical potential and isospin chemical potential with the chemical potentials chosen so the that 
baryon density is twice the isospin density. It is unclear whether the x?T techniques used here can be generalized 
easily to deal with this problem since they are based on the assumption that the only light scale in the problem is 
the pion mass. On the other hand, the inverse of the nucleon-nucleon scattering length, I/a, defines a mass scale 
much lighter than to^ || . This suggests that one should develop a systematic expansion treating both 1/a and is 
light. Formally, such expansions have been developed for the two nucleon problem Q however, in that context there 
is clear evidence for problems with convergence of the chiral expansion JTo| ] . However, the effects of chiral symmetry 
may be very different in the case of neutron matter from those in the N-N scattering case and a viable expansion may 
prove possible. The possibility that at some density neutron matter exhibits pion condensation is interesting since 
the fact that the propagation of nucleons through a pion condensate changes characteristically above a certain isospin 
chemical potential could play a significant role. 

With the above caveats, we have looked at the mean-field Fermi gas in the present framework. The possibility of 
pion condensation arises from the term — i/i?/^ sin 2 a in the energy, which for sufficiently large p z can overcome the 
mass term —to 2 / 2 cos a; the question is whether a gas of fermions tends to stabilize or destabilize the phase with 
condensed pions. We see from eq. ( |l3| ) that for a given value of p l7 the lowest single-nucleon energy is at a — 0, the 
normal vacuum. Analytic calculation shows that for fixed p 1 and baryon density, a = is also a minimum of the full 
fermion energy, whereas a — ^ir is a maximum. Thus the fermions do not favor the pion-condensed phase. Numerical 
calculations were done for the case with the isospin density constrained to equal minus half the baryon density, as 
in neutron matter. In the chiral limit a = ^ir (which minimizes the pion energy) is a minimum of the full pion and 
fermion energy only for baryon densities below about a third of nuclear matter density (po)- For higher densities the 
minimum moves round the chiral circle, reaching a = at about two-thirds of po. Thus even in the chiral limit the 
condensed phase is a low-density phenomenon. With the physical pion mass, pion condensation requires p 1 > m v at 
zero baryon density; in neutron matter such a value is reached at around 4po- Thus the Fermi gas model does not 
predict this form of (s-wave) pion condensation, much less the new phase discussed in this paper, in neutron stars. 
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